در این مقاله چهارچوب تئوریک یک مدل شکست ترموهیدرومکانیکی (مدل شکست (THHMD برای محیط متخلخل غیراشباع غیرهم دما ارائه شده است. چهارچوب این مدل بر اساس استفاده از متغیرهای حالت مستقل (تنش خالص، مکش و تنش حرارتی) می باشد. قانون رفتاری شکست از مفاهیم میکرومکانیکی و توضیح پدیده ای نشات گرفته است. روابط مزدوج ترمودینامیکی تنش- کرنش از انرژی آزاد که جمع انرژی های تغییرشکل الاستیک شکست و پتانسیل های کرنش باقی مانده است گرفته شده است. سختی های مکانیکی شکست با اعمال اصل انرژی الاستیک معادل (PEEE)برای هر هر متغیر حالت تنش محاسبه می شود. اثر شکست بر انتقالات آب و بخار با معرفی نفوذپذیری ذاتی ناشی از شکست اعمال شده است. مدل THHMD شکست در برنامه المان محدود q-Stock وارد شده است.در این مقاله پس از توضیحات مدل مذکور یک مطالعه عددی برای بررسی اثرات ناحیه شکسته ناشی از حفاری (EDZ) به پاسخ مدفن زباله های اتمی در طول فاز حرارت دهی صورت گرفته است. روند مشاهده شده با انتظارات تئوریک تطابق مناسبی دارد.